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Abstract. We investigate here the Cooper pairing of fermionic atoms with mismatched fermi surfaces using 
a variational construct for the ground state. We determine the state for different values of the mismatch of 
chemical potential for weak as well as strong coupling regimes including the BCS BEC cross over region. We 
consider Cooper pairing with both zero and finite net momentum. Within the variational approximation 
for the ground state and comparing the thermodynamic potentials, we show that (i) the LOFF phase is 
stable in the weak coupling regime, (ii) the LOFF window is maximum on the BEC side near the Feshbach 
resonance and (iii) the existence of stable gapless states with a single fermi surface for negative average 
chemical potential on the BEC side of the Feshbach resonance. 
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1 Introduction 

Fermionic superfluidity driven by s-wave short range in- 
teraction with mismatched fermi surfaces has attracted 
attention recently both theoretically [1] as well as experi- 
mentally [2]. The experimental techniques developed dur- 
ing the last few years, at present, enable one to prepare 
atomic systems of different compositions, densities, cou- 
pling strengths as well as the sign of the coupling in a 
controlled manner with trapped cold fermionic atoms us- 
ing techniques of Feshbach resonance [3]. Because of the 
great flexibility of such systems, the knowledge acquired 
in such studies shows promise to shed light on topics out- 

, side the realm of atomic physics including dense quark 
matter in the interior of neutron stars [4], high Tc super- 

, conductors [5] as well as for the physics of dilute neutron 

■ gas [6]. 

\ In degenerate fermi systems, it is well known that an 

■ attractive interaction destabilizes the fermi surface. Such 
' an instability is cured by the BCS mechanism character- 
ized by pairing between fermions with opposite spins and 
opposite momenta near the fermi surface leading to a gap 
in the spectrum. In this system the number densities and 
the chemical potentials of the two condensing fermions 
are the same. On the other hand, there are situations 
where the fermi momenta of the two species that condense 
need not be the same. Such mismatch in fermi surfaces 
can be realized in different physical systems e.g. super- 
conductor in an external magnetic fleld or a strong spin 
exchange field [7,8], mixture of two species of fermionic 
atoms with different densities and/or masses [2] or charge 
neutral quark matter that could be there in the interior 



of neutron stars. The ground state of such systems shows 
possibilities of different phases which include the gapless 
superfluid phase [9-11] with the order parameter as non 
zero but the excitation energy becoming zero, the LOFF 
phase where the order parameter acquires a spatial varia- 
tion [7,12-14]. The LOFF phase might also manifest in a 
crystalline structure [15]. 

When the mismatch in chemical potential is small as 
compared to the superfluid gap, the superfluid state does 
not support the gapless single particle excitations. When 
the chemical potential difference is large as compared to 
the superfluid gap there could be possibility of breached 
pairing superfluidity [9,10], where there are two fermi sur- 
faces with the single particle excitation vanishing at two 
values of the momentum. When the mass difference be- 
tween the two species is large it could also lead to interior 
gap superffuidity where the quasi particle excitation van- 
ishes at a single momentum [17]. The stability of such 
configurations has been studied demanding positivity of 
the Meissner masses [14, 16] as well as the number suscep- 
tibility [18]. Such an analysis has also been done for rela- 
tivistic fermions with four fermi interactions [16]. In these 
analyses the ground state was not considered explicitly. In- 
stead, the chemical potential and the gap functions were 
treated as parameters, with the whole parameter space 
being analyzed systematically for stability criteria which 
were related to different possible quasi particle dispersions 
[16,18]. 

Our approach to such problems has been variational 
with an explicit construct for the ground state [10,19-22]. 
The minimization of the thermodynamic potential with 
respect to the functions used to describe the ground state 
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decides which phase would be favored at what density 
and/or coupUng. This has been apphed to quark matter 
with charge neutrahty condition [21-23] as weh as to the 
systems of cold fcrmionic atoms [10]. In the present paper 
we would like to extend the method to have the ansatz 
state general enough to include the LOFF state with a 
fixed momentum for the condensate. Thus our approach 
will be complementary to that of Ref. [18] in the sense that 
we shall use explicit state and solve the gap equation to de- 
termine gap function. Comparison of thermodynamic po- 
tential for different phases after the gap equation is solved 
will decide which phase will be favored at what coupling, 
density and the mismatch in fermi momentum. We also 
take the analysis further as compared to that in Ref. [18] 
to include the LOFF state, with a single plane wave. 

For the case of vanishing mismatch in chemical poten- 
tials of the pairing species, a functional integral formula- 
tion with a saddle point approximation was considered to 
describe the crossover from BCS to Bose superfluidity [24, 
25] . A systematic analysis of the effect of gaussian fluctua- 
tions around the saddle point solution was also performed. 
It was shown that while the corrections from the fluctu- 
ations are extremely important for temperatures close to 
the critical temperature, particularly for strong coupling, 
the same is very small for all couplings for temperatures 
small compared to the critical temperature. The ability of 
this saddle point solution which at first sight might have 
been expected to work only for small couplings, to repro- 
duce the strong coupling Bose limit is reassuring. This 
gives some confidence in its validity for intermediate cou- 
pling results where no obvious small expansion parameter 
is known [24]. As we shall show in the following, the re- 
sults arising from our ansatz for the symmetric case of 
equal chemical potentials corresponds to this saddle point 
solution. 

We organize the paper as follows. In section 2, we 
briefly introduce the model and the discuss the ansatz for 
the ground state. In subsection 2.1, we compute the expec- 
tation value of the Hamiltonian with respect to the ansatz 
state to calculate the thermodynamic potential. Minimiz- 
ing the thermodynamic potential, we determine the ansatz 
functions in the ground state and the resulting superfluid 
gap equations. In section 3, we first discuss various pos- 
sible phases and the corresponding results of the present 
investigation. Finally, we summarize and conclude in sec- 
tion 4. 



where, ip°-{r) is the annihilation operator of the fermion 
species 'a', and, g is the strength of interaction between 
the two species of the fermionic atoms. Throughout the 
paper wc shall work in the units of ?i = c = fcg = 1. 

While considering Cooper pairing between fermions 
with different fermi momenta, it was realized by Larkin 
- Ovchinnikov and Fulde- Ferrel (LOFF), that it could be 
favorable to have Cooper pairing with nonzero total mo- 
mentum unlike the usual BCS pairing of fermions with 
equal and opposite momentum. We shall examine here 
whether an ansatz with a condensate of Cooper pair of 
fermions having momenta k + q/2 and — k + q/2, thus 
with a nonzero net momentum q of each pair is favored 
over either BCS condensate or the normal state without 
condensates. Here, our notation is such that k specifies a 
particular Cooper pair while q is a fixed vector, the same 
for all pairs which characterizes a given LOFF state. Here 
the magnitude of the momentum |q| shall be determined 
by minimizing the thermodynamic potential while the di- 
rection is chosen spontaneously. 

With this in mind, to consider the ground state for this 
system of fermionic atoms with mismatched fermi surfaces 
(/ii ^ 112), we take the following ansatz for the ground 
state as 

\n)=Ug\vac) (2) 

where, \vac) is annihilated by V'r(r) and the unitary op- 
erator Uq is given as 

Ug = exp (^J r (k + |)V''(-k + |)te„6/(k)dk - h.c.^ . 

(3) 

In the above, k identifies a fermionic pair and each pair 
in the condensate has the same net momentum q. The 
function /(k) is a variational function to be determined 
through the minimization of the thermodynamic poten- 
tial. In the limit of zero net momentum (q = 0) this ansatz 
reduces to the BCS ansatz considered earlier [10,22,23]. 
Such an ansatz breaks translational and rotational invari- 
ance. In position space, as we shall show that it describes 
a condensate which varies as a plane wave with wave vec- 
tor q. To include the effect of temperature and density 
we write down the state at finite temperature and den- 
sity, ]/2(/3,/Lt)), taking a thermal Bogoliubov transforma- 
tion over the state | i?) using thermofield dynamics (TFD) 
method as [26,27] 



2 The ansatz for the ground state 

To discuss the superfluidity for fcrmionic atoms, wc con- 
sider a system of two species of non-relativistic fermions 
with chemical potentials /xi and 1J.2 and having equal masses. 
These two species can e.g. be the two hyperfine states of 
■^^Ca. Further we shall assume a point like interaction be- 
tween the two species. The Hamiltonian density is given 
as 

W(r) = V>„(r)t(-|^-/,„)V.„(r)+ff^t(r)Vi(r)V2(r)V2(r) 

(1) 



where U^^n is given as 



with 



Bt(/3,^)= J [/"(k)t^«(k,/3,M)/''(-k)t 



dk. 



(4) 



(5) 



(6) 



In Eq.(6) the ansatz functions 6*" (k, /3, /i), as we shall see 
later, are related to the distribution function for the a-th 
species of fermions, and, the underlined operators are the 
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operators in the extended Hilbert space associated with 
thermal doubling in TFD method. All the functions in 
the ansatz in Eq.(4), the condensate function /(k), the 
distribution fimctions 0*(k, /?, /i) are obtained by the min- 
imization of the thermodynamic potential. We shall carry 
out this minimization in the next subsection. 



2.1 Evaluation of thermodynamic potential and the 
gap equation 

Having presented the trial state Eq.(4), we now proceed 
to minimize the expectation value of the thermodynamic 
potential with respect to the variational fimctions /(k), 
0"(k, /?) in the ansatz. Using the fact that, the variational 
ansatz state in Eq.(4) arises from successive Bogoliubov 
transformations, one can calculate the expectation values 
of various operators. In particular, for the condensate we 
have 



{^2{|3,^Ji)\r{v)i,\v)\Q{|3,^l)) 

e"*" exp(iq • r) 



sin 2 



2(27r)3 

X [sin^ e^{-\<i + q) - cos^ r (k)] dk 
= e"''cxp(z:qT)Jg, 



(7) 



which for nonzero Iq describes a plane wave with a wave 
vector q. Once one has demonstrated the instability to 
formation of a single plane wave, it is natural to expect 
that the state which actually develops has a crystalline 
structure. LOFF, in fact, argued the favored configura- 
tion to be a crystalline condensate which varies in space 
like a one dimensional standing wave cos(q • r) with the 
condensate vanishing at the nodal planes. Which crystal 
structure will be free energetically most favored is still to 
be resolved. The present variational ansatz Eq.(4) is a first 
step in this direction to decide which crystalline structure 
will finally be chosen through minimization of the ther- 
modynamic potential. 

Similarly the number operator expectation value for 
the atoms of species 'a' is given as 



p« = M/3,M)ir(r)V(r)|r2(/3,A^)) 
= [cosV(k-f)sin2r(k) 

+ le"*-] sin^ /(k - |) cos' ^*'(-k + q)] dk. 



(8) 



which is independent of the space coordinate r. 

Using Eq.s (7) and (8) it is straightforward to calculate 
the expectation value of the Hamiltonian given in Eq.(l). 
We obtain 



e={n{|3,^i)\H\Q{|3,^l) 



cosV(k)sin2 0i(k+5) 



sinV(k) cos' 02(-k+^) 



+ (e(-k+|)-M2) 



cos' /(k) sin' 02(-k + ^) sin' /(k) cos' ^i(k + |) 
+9II + 9P1P2 



(9) 



where, the first term arises from the expectation value 
of the kinetic part of the Hamiltonian and the last two 
terms arise from the four point contact interaction term. 
Here, e(k) = k'/(2m), and, the quantities Iq, p° (a = 
1,2) are defined in Eq. (7) and Eq. (8) respectively. The 
thermodynamic potential is given as 



f2 : 



1 



(10) 



where, we have denoted s as the total entropy density of 
the two species given as [26,28] 



E 



dk (sin' 6*" ln(sin' 6*") + cos' O'' ln(cos' 0"" 

(11) 

We now apply the variational method to determine 
the condensate function /(k) in the ansatz Eq.(4), by re- 
quiring that the thermodynamic potential is minimized : 
df2/df{k) = 0. Such a functional minimization leads to 



tan2/(k) = 



9lq _ ^ 



Co - V 



(12) 



In the above, = (e(k + f ) + e(-k + a))/2 = (k' + 

q'/4)/(2m) is the average kinetic energy of the two con- 
densing fermions. Similarly, 9 — (yi +z^2)/2, is the average 
of the interacting chemical potential of the two fermions 
with v"- = jji"- — gp'^ and 67 = Eg — v. Further, we have 
defined A = —glq, with, Iq as defined in Eq.(7). It is 
thus seen that the condensate function depends upon the 
average energy and the average chemical potential of the 
fermions that condense. Substituting the solutions for the 
condensate functions given in Eq. (12) in the expression 
for Iq in Eq.(7), we have the superfluid gap equation given 

by 



A = — 



(27r)3 



dk- 



2a;„ 



1-sin' 6ii(k+^) 



sin'^(-k+|) 



(13) 



In the above, 



g — Y + and the thermal functions 
0°(k, /?) are still to be determined. 

The minimization of the thermodynamic potential with 
respect to the thermal functions ^a(k) gives 



sin' 



exp(/?a;a) -h 1 



(14) 



with the quasi particle energies given as a;i(k-|- ^) = Wg -|- 
and t<;2(-k + |) = a)g - Jj, with (5^ = (^1 - ^2)/2 = 
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(e(k + f ) - e(-k + a))/2 - (ly^ - z.2)/2 = - 5,. It is 
clear from the dispersion relations that it is possible to 
have zero modes, i.e., w° = depending upon the values 
of 5f and 5,j. So, although wo shall have nonzero order 
parameter A, there can be fermionic zero modes, the so 
called gapless superconducting phase [29,30]. 

Now using these dispersion relations, and the super- 
fluid gap equation (Eq.(13)), the thermodynamic poten- 
tial (Eq.(lO)) becomes 



zero temperature and zero density contribution. The reg- 
ularized gap equation is then given as 



(2 



^^ln(l + exp(-/3a;„)) 



Z\2 

— (15) 



Q = 



+ 



1 



(2' 



^ j dk[a;i0(-a;i)+a;20(-a;2)]-y.(16) 



To compare the thermodynamic potential with respect 
to the normal matter, wc need to subtract out the zero 
gap and zero momentum (q = 0) part from the thermo- 
dynamic potential given in Eq. (16). This is given as 

Q = I2(q, Zi, 5) - i7(q = 0, zi = 0, (5) 



+ 



(27r) 



^ j rfk[a;i0(-a;i)-Fa;20(-a;2)] 



(27r) 
1 



(17) 



where, we have denoted ^ = k^/(2m) — the excitation 
energy with respect to the average fermi energy for normal 
matter. Further, we have assumed, without loss of gener- 
ality, 5,^ > 0. The supcrfluid gap in Eq.(17) satisfies the 
gap equation Eq.(13) which at zero temperature reduces 
to 



1 _ 



^ f dk-^ i-e( 



■a;i(k+|) 



-k + 



(18) 



The above equation is ultraviolet divergent. The origin of 
this divergence lies in the point like four fermion interac- 
tion which needs to be regularized. Wc tackle this problem 
by defining the regularized coupling in terms of the s-wave 
scattering length as was done in Ref. [24,31] to access 
the strong coupling regime [10,33] by subtracting out the 



1 



M 



9r 47ras 



(27r)3 / 



dk 



1 - e{-uji)0{-uj'2) 

2u)a 



2e(fc) 



(19) 



Using Eq.(18)and Eq.(19) in Eq.(17) to eliminate g in 
favor of qr, one can obtain 



1 



In what follows we shall concentrate on the case for zero 
temperature. In that case the distribution functions, Eq.(14) 
become - functions i.e. sin^ = 0{—LOa)- Further, us- 
ing the identity limo_>oo lii(l -|- exp(— ax)/a = —x0{—x), 
in Eq. (15), the zero temperature thermodynamic poten- 
tial becomes 



+ 



(27r)3 
1 



(27r)3 



Z\2 Z\2 

(Wl - 7Z^)0{-OJl) + - — )6)(-W2) 



2cj„ 



\\ + k)Q{-\U-k) 



+ (l^",l-^^)e(-l^,l + <5«))dk. 



(20) 



where, |^q| = (k^ -|-q^/4)/2m- and 8^ = k- q/2m- (5;,. 
We might note here that the thermodynamic potential 
difference as in Eq.(20) is cut-off independent and is fi- 
nite. The only other quantities needed to calculate the 
thermodynamic potential in Eq.(20) arc the chemical po- 
tentials: the average chemical potential of the two species 
u = {vx+v2)/2 and their difference 5^ = {vi — ^2)- These 
two quantities can be fixed from the average number den- 
sities as 

p=\{{4'4>i) + {4i^2)) 

and the difference in number densities 

5p = ^((V'l^i) - (V'IV'2)) 
1 



2^ 

1 

2(2^ 



j d\<i{0{-UJi)-0{-UJ2))^ 



(22) 



which does not depend on the condensate function explic- 
itly. In the following we shall discuss which phase is ther- 
modynamically favorable at what density as the chemical 
potential difference is varied for a given coupling and a 

given average density. 

For numerical calculations, it is convenient to write the 
Eq.s (19-22) in terms of dimcnsionlcss quantities. Thus we 
make the substitutions |k| = kfX, q = fc/y, A = CfZ, V = 
e/v, = EfS, where, kf is the average fermi momentum 
and € f is the corresponding fermi energy. In terms of these 
dimensionless variables, the gap equation Eq.(19) reduces 
to 



Snkfa 



1 

4^ 



dxx 



v/(x2 + 2/2/4 - vf + 
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1 

1 



dxx^ 



dt{0{-UJi) + e(-Cu2)). 



(23) 



Here, the quasi particle energies (2)a = 1:2) in units of 
average fermi energy are 

lDi — ujy + xyt ~ 5, UJ2 — ojy ~ xyt + 5 (24) 



where, Cjy — ^ (x^ + y2/4 — v)"^ + ^ where, v is the aver- 
age of the excitation energies of the two quasi particles in 
units of fermi energy. Further, in Eq.(23), the integration 
variable t = cos 9 and, we have chosen the momentum q 
of the condensate to be in the z- direction. 

Similarly, the average density equation Eq.(21), in units 
of Kj I dn^ , reduces to 



1 



dxx^ 



dt 



(1 



(25) 



Finally, the thermodynamic potential difference be- 
tween the condensed phase and the normal phase can be 
expressed in terms of these dimensionless variables as 




27r2 



X dx- 



dt 



x^dx 



y_ 

4 



2t2), 



Z 

20}. 

2d), 



■)e(-'^2) 



Fig. 1. Superfluid gap as a function of coupling 1/kfa. The 
expoential rise of the gap in the BCS side {kfa < 0) is as 
expected from weak coupling results. 



and, the number density equation, Eq.(25), becomes 



-- J x'dxdtl^i\Q+S^^ 



3 

2./0 



dxx 1 



(28) 



e/4 
27r2 



Lyl 



ly2 



iy^ ^ (26) where, z and v are respectively the superfluid gap and 

average chemical potential in units of fermi energy. 
5 and Iy^{i = 1, 2, 3) are the three in- K is easy to show, via integration by parts, that the 



where, S^^^ = xyt 

tegrals in the above equation.^ Eq.s (23)', (25) and Eq.(26) integrals h and Ip can be rewritten as 
will form the basis for our numerical investigation to de- 
cide regarding the phase structure. Before going to the 
numerical evaluations let us discuss the different phases. 



1 



= ^ z) - z'^J{v,z)) 



and 



3 Results 

3.1 Homogeneous symmetric phase 



Ip = z'^{I{iy,z) + i^J{iy,z)) 



(29) 



(30) 



where, we have defined the integrals /(t^, z) and J(i^, z) as 



Let us first discuss the symmetric homogeneous phase 
corresponding to vanishing chemical potential difference 
{di, = 0). In that case the renormalized gap equation, ^^^^ 
Eq.(23), reduces to the usual BCS gap equation [10,24] 



J{iy, z) 



dx 
uj{x) 



X dx 
uj{x)^ 



(31) 



(32) 



1 



1 

4^ 



dxx 



^ [xP- — I/)2 + 2;2 



I with ijj{x) ~ [xp- — 1^)2 + z2 

/ The nice thing about these integrals l(y^ z) and J(z^, z) 

(27) is that they can be expressed in terms of elliptic integrals 
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Fig. 2. Self consistent solution for the chemical potential in 
units of fermi energy as a function of coupling l/kfU. 



[25]. They can be written as 



This leads to the weak couphng BCS limit for the gap as 
z = 8e~^ exp(— 7r/(|fcja|)) in units of fermi energy. Simi- 
larly in the BEC limit (l/(fc/a) oo), one expects tightly 
bound pairs with binding energy Eb = l/(ma^) and non- 
degenerate fermions with a large negative chemical po- 



tential [24]. In this limit 



2/(41/-^) << 1. Taking 



this asymptotic limit of the elliptic integrals, the chem- 
ical potential v is half of the pair binding energy. The 
corresponding gap is z = (16/7r)^/^/-^/fc/a . 

As may be seen from Fig. 2, the chemical potential 
changes sign at Kf a ^ 0.56 signaling the onset of the BEC 
phase. Near to the unitary limit e.g. for fc/a = 0.01, the 
ratio of condensate to the chemical potential turns out 
to be A/h' ~ 1.195, while the chemical potential itself is 
numerically evaluated to be J^/e/ = 0.58. While the ratio 
of gap to the chemical potential A/i/ is very close to the 
value A/v = 1.2 obtained through quantum Monte Carlo 
simulations, the value of the ratio of chemical potential 
to the fermi energy, ly/ef, is higher as compared to the 
value I'/ef — 0.42 obtained from quantum Monte Carlo 
simulation [34]. Such results arising from the variational 
calculations may be compared with the corresponding val- 
ues obtained from other nonperturbative calculations like 
e- expansion [36] which turn out to be A/v ~ 1.31 and 
I'/cf = 0.475. We might note here that, a systematic anal- 
ysis of the gaussian fluctuation around the saddle point ap- 
proximation was considered in Ref. [24]. The corrections 
to the saddle point approximation were seen to be small 
for temperatures small compared to the critical tempera- 
ture, for all couphngs [24]. 



and 



J{v,z) 



2(^.2 + ^2)1/4 
(z2+j.2)l/4 



F{tt/2,k) 



(33) 



F{TT/2,k) 



2(z2 + iy2)l/4(^ + yj2-^) 



(34) 



where, F{^,k) and E{^,k) are the elliptic integrals of 
first kind and second kind respectively. In the above, 

= {i/ + Vz2+T^)/(2vzM-i^). The thermodynamic 
potential difference Eq.(26) reduces to 



n{6 = 0,z) = 



27r2 



dxx 



X — V\ — LO - 



z 

2(b 



j)(35) 



where, Cj — ^ (xp- — jy)^ + z^. 

In figures 1 and 2 we plot the superfluid gap and the 
chemical potential as functions of coupling 1/fcya obtained 
by solving the coupled gap equation (Eq.(27)) and num- 
ber density equations (Eq.(28)). In the weak couphng BCS 
limit (l/(A:/a) —oo), we have v » z and the modu- 
lus parameter ~ 1 — z"^ j {'\,v'^') is close to unity in the 
argument of the elliptic integrals. Expanding the elliptic 
integrals in this limit, the gap equation Eq.(29) reduces to 



1 



47r2 



In 



8i^ 



(36) 



3.2 Isotropic superfluid with mismatch in chemical 
potentials 

With a nonzero mismatch between chemical potentials i.e. 
Si, 7^ 0, let us consider superfluidity with a homogeneous 
condensate i.e. condensate with momentum q = 0. In that 
case the dispersion relations (in dimensionless variables) 
of the quasi particles (Eq.(24)) are given as 



^(a;2 - iy)2 + z2 _|_ ^. i2)2 = \/{x'^ - vY + z2 



(37) 

where, v, z and 5 are respectively the average chemical 
potential, the gap and the difference in chemical potentials 
in units of fermi energy. 

For any value of the (dimensionless) gap 2; > 0, the 
average chemical potential u and the difference in chem- 
ical potentials 5 > Q, the excitation branch loi has no 
zeros similar to the case of a usual superconductor. The 
quasi particles at the fermi surface have finite excitation 
energy given by the gap that is enhanced here by the mis- 
match 5. On the other hand, the excitation branch uj2 
can become zero depending upon the value of the gap z, 
the average chemical potential u and the mismatch pa- 
rameter 5. Indeed, solving for the zeros for uj2 one ob- 
tains that UJ2 vanishes at momenta (in units of fermi mo- 
mentum), a;2 , ^ u ± ^5"^ - z2. Clearly, this has 



c2 , . ^ u±^/6' 

max/ mm ^ 

no solution for 5 < z. Further, for the case oi 5 > z, 
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there will be no solutio ns if the average chemical poten- 
tial V is smaller than —\/5'^ — since both a;^^^. and a;^j„ 
will become negative. On the other hand, if 5 > z and 
— \J5'^ — < D < \/S^ — z^, then the expression for x\^^^ 
is positive while that of a;^j„ is negative. Thus there is 
only one zero for a;2. This case will correspond to 'interior 
gap' solution, where, the unpaired fermions of the first 
species occupy the entire effective fermi sphere. Finally, 
when x^Q^ > and a;^j„ > 0, there are two zeros for uji. 
This will correspond to breached pairing case where un- 
paired fermions of first species occupy the states between 
the two fermi spheres decided by Xmax and Xmin- 

With nonzero mismatch in the chemical potential, the 
gap equation for the homogeneous phase is given as 



3.0 



Snkja 



1 
1 

Iz-Is 



dxx 



1 



1 



L0{x) 



uj{x) 



dx 



Here, we have introduced the notation ui{x) —\/ (a;^ — vY' + 
as the average of the quasi particle energies of the two 
species with v as the average chemical potential in units 
of fermi energy. 

Further, the thermodynamic potential of Eq. (26) re- 
duces for homogeneous condensates to. 




n 



27r2 

+ / x^dx ( (oj 



X dx I \x~ 

z2 



uj(x) 



z 

2Q 



2uj{x) 
x^dx0{S~ \x^ ~ v\ 



S). 



~ 27r2 



(39) 



where, we have defined the three integrals appearing in 
Eq.(39) as /i, I2 and I3 respectively. The integral I3 can 
be evaluated directly as 

h = J x^dx0{6 - \x^ - v\){\x^ -v\-^) 

= ^/^^ ~ ^(e(^+)^r + e(z._)^^/^) (40) 

with v± — V :Lb. 

Similarly, the integral I2 can be rewritten as 



I2 = J x'^dx{uJi 



2Q{x) 



x^dx I uj{x) — S - 



2u{x) 



(41) 



-1.0 -0.5 0.0 
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Fig. 3. Ratio of critical chemical potential difference to the 
gap as a function of coupling strength 1/kfa. Gapless phase 
appears for coupling 1/kfa > 1. The region between the dotted 
and the dashed line shows the gapless phase. 



is characterized by the average chemical potential being 
equal to the fermi energy so that 1^=1 and is much 
larger than the magnitude of both the chemical poten- 
tial mismatch S in dimensionless units and the gap z, also 
in dimensionless units. Thus this weak coupling analysis 
excludes the scenario with single fermi surface i.e. with 
the excitation energy 102 having a single zero. Further, 
we here compare the thermodynamical potential for fixed 
chemical potentials. It is expected that the same analysis 
with fixed number densities will lead to different results 
- namely a mixed phase scenario which is an inhomo- 
geneous phase where a fraction of the space is in normal 
phase while the remaining is in the BCS phase [33]. In this 
weak coupling limit, we can have ordinary superconduc- 
tivity without gapless modes or can have breached pairing 
with LU2 becoming zero at two values of the momentum. 
The latter situation will arise for much larger values of the 
mismatch parameter S as compared to the gap z. In the 
weak coupling limit, the integral Ii becomes 



h - 



v^dx 



- uj{x) 



z 

2Q 



(42) 



Before going to the numerical solution for the homo- 
geneous condensates with mismatch in chemical poten- 
tial, let us discuss the weak coupling BCS case. This limit 



'v = x„iax, SO that, I2 — 0. Fi- 



In the same limit, Xmin- 
nally, the integral ~ ^JvS^ /A. Thus in the weak coupling 
limit, the thermodynamic potential difference between the 
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condensed and the normal phase is 

r2^^(2<52_z2) (43) 

For the ratio of the chemical potential difference to the 
superfluid gap 5/z < 1/a/2, the Clogston-Chandrasekhar 
limit [32] , the thermodynamic potential difference becomes 
negative. In this regime of 5/z, the gap is the same as the 
(5 = gap and the density remains the same even though 
5 is nonzero. This critical chemical potential difference be- 
low which BCS pairing is possible depnds on the coupling 
strength. 

In Fig. 3 we plot, as a function of the coupling 1/ Kfa, 
the quantity (5c/^ the ratio of the maximum chemical po- 
tential difference to the gap that can sustain pairing. For 
the region below the solid and the dotted line corresponds 
to the parameter region where BCS pairing is the state of 
minimum thermodynamic potential. As may be seen from 
the figure, for weak coupling, the critical chmecal poten- 
tial difference (in units of superfluid gap) approaches the 
Clogston-Chandrasekhar limit. As coupling is increased 
from BCS to BEC side 5c/ A increases monotonically. 

We do not flnd any breached pairing solutions in the 
BCS region with lower thermodynamic potential as com- 
pared to the normal matter. However, we do find the in- 
terior gap solution in a range of 5^ in the strong cou- 
pling BEC regime with negative average chemical poten- 
tial. This happens for l/kfa > 1. This is given by the 
region between the dotted and the dashed line in Fig. 3. In 
this region the density difference between the two species 
is also nonzero. In the region of the parameter space that 
is above the solid line or to the left of the dashed line is the 
region where no pairing is possible and the normal matter 
is the state of minimum thermodynamic potential. 

In Fig. 4 we plot the gap as a function of difference 
in chemical potential for l/kfa = 2 as a typical value 
for strong coupling BEC regime where gapless phase ex- 
ists. In this region, the thermodynamic potential differ- 
ence between the condensed phase and the normal phase 
is negative. The value of the gap decreases continuously 
from the symmetric BCS phase to the gapless phase as 5^ 
is increased leading to nonzero difference in the number 
densites. The transition from gapless phase to the normal 
phase which occurs at 5i, ~ 3.85Z\o is however discontin- 
uous. In the gapless phase the density difference between 
the two species is non vanishing. 

As mentioned earlier, in our calculations we do not 
keep the difference in the densities of the two species fixed. 
The calculations are performed with a fixed average den- 
sity and given value of chemical potential difference. In the 
gapless phase, the density difference comes out to be non- 
vanishing. In Fig. 5 the dependence of the gap on the den- 
sity difference 5p is shown for the BEC region where gap- 
less phase exists. We have taken here the value of I/k/o 
as equal to 2. Superconductivity is supported for a max- 
imum density difference of 5p ~ 0.52 in units of k^/Gir'^ 
beyond which the system goes over to the normal matter 
with zero gap. For coupling l/fcja<l,wedo not find any 
superfiuid phase free energetically favorable for any non 



2.0 




0.5 



0.0 ' ■ ■ ■ • ■ ■ • ^ 

12 3 4 

6/Ao 



Fig. 4. Superfluid gap as a function of the difference in the 
chemical potentials of the two species. The curve is plotted for 
l/kfa = 2. 

zero value for the density difference, although a chemical 
potential difference can still support a Cooper paired BCS 
phase with zero net momentum.. 

In Fig. 6, we have plotted the phase diagram in the 
plane of average chemical potential and the chemical po- 
tential difference, both the quantities being normalized 
with respect to the superfluid gap. The region above and 
to the right of the solid line correspond to positive ther- 
modynamic potential and is not stable. The region below 
and to the left of the solid line shows Cooper pairing. The 
region between the dashed and the solid line shows the 
gapless phase. The vertical line at 5/ A = l/\/2 = 0.717 
indicates the Clogston Chandrasekhar limit. These results, 
corresponding to explicit solutions for the gap as a func- 
tion of coupling are in conformity with the results obtained 
in Ref. [18]. 



3.3 The LOFF Phase: anisotropic superfluid with q 7^ 

When the value of the difference in chemical potentials, 5, 
exceeds the Clogston-Chandrasekhar limit, it can be free 
energetically favorable to have condensates with nonzero 
momentum q, given as in Eq.(7). The dispersion relations 
for the quasi particles are given as 

ti;i(k) = cZ'(k) -|- 5^ 
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Fig. 5. Superflmd gap as a function of difference in number 
densities of the condensing species. The curve is plotted for 
l/kfu = 2. 



2m 



(44) 



W2(k) = w(k) - 



+ 



kg 

2m 



k2 + :i- ) /2m - 



(45) 



The corresponding gap equation, number density equation 
and the thermodynamic potential are already given in Eq.s 
(23), (25) and (26) respectively. 

As in the previous subsection, before going to discuss 
the detailed numerical solutions for strong couplings, let 
us consider the weak coupling limits of the above equa- 
tions which can give an insight to the numerical solutions 
that can be obtained in the appropriate limit. In this limit, 
the average chemical potential is same as the fermi energy 
sothaty^ = I'/ef ^ 1 and is much larger than the mo- 
mentum (in units of fcrmi momentum) y, gap ( in units 
of fermi energy) z and the asymmetry in chemical poten- 
tial (in units of fermi energy) , d. Further, the excitation 
energies in units of fermi energy can be approximated as 



^1.2 — 



= UJ,, ± St 



(46) 




6/A 



Fig. 6. The ratio of average chemical potential to the gap as a 
function of difference of chemical potential The region between 
the dotted and the solid line shows the gapless phase. 



where, ujy ~ y + Av{x — y^)^, and, St — yt — S = 

5{pt — 1), with p — y/5. Just to remind ourselves, y is 
the condensate momentum in units oi kf and S is half of 
the difference in the chemical potentials of the two species 
in units of fermi energy £/ = k'j/2m. Using Eq.(36), one 
can eliminate the coupling 1/kfa from the gap Eq.(23) to 
obtain the weak coupling LOFF gap equation as 



In 



(?) 



dxdt 



iOi-LOl) + 0(-o.2)) . (47) 



The theta functions in the integral above restrict the limits 
of X integration to be in the range (a;,„i„, a;m(ja;),with. 



1 



2V^ 



2U/2 



(48) 



Performing the a:- integration within the limit, the gap 
equation, Eq.(47), reduces to 



ln(f^^ 



= - / dte{\St\ 



1 In 



\St\ + ^W^z 



\st\ - ,/sr^z 



(49) 



Again, the theta function restricts the angular integration 
where either or both the quasi particle excitation energies 
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(aji^2) become gapless. This occurs for the case when the 
gap z is less than (1 + as otherwise the domain of inte- 
gration specified by the theta function 6>(|<5t| —z) becomes 
null. Depending upon the value of the gap z as compared 
to (p — 1), either one of the two quasi particles or both 
become gapless [37]. Using this fact, the integration over t 
can be performed in Eq. (49) leading to the weak coupling 
LOFF gap equation as 



In 



\ z 



1 



4p 



In 



1 — xi 
1 + 



2x1 



p-1 
4p 



In 



= 0, 
where, xi, 



xi are parameters 



xi = 6) 1 



and. 



2^2 



1 



(P+l)'5 



\p-\\8 



1 — a;2 
l + a;2 



1/2 



2x' 



1 - 



(52(p+l) 



1 



<52(p-l)2 



1/2 



(50) 



(51) 



(52) 



We note here that the gap equation Eq.(50), is identical to 
the one derived in Ref.[8] and Ref.[37] for LOFF phase in 
the superconducting alloy with paramagnetic impurities 
and in quark matter respectively. 

Next, let us look at the thermodynamic potential Eq.(26) 
in the weak coupling limit. The thermodynamic potential 
in Eq.(26) can be written as 



n ■ 



27r2 



Ly2 



(53) 



where, /yi(i=l,2,3) are the three integrals of Eq.(26). It 
is easy to show, using Eq.(42), that in the weak coupling 

limit 

t/V4)z2 



(54) 



The evaluation of the integral Iy2^ is similar to the evalu- 
ation of the integral on RHS of the gap equation Eq.(47). 
The result is 



ly2 



12 



[p-lf 



(55) 



P P 

The integral ly^ is straightforward to evaluate and is given 
as 




yl8 



Fig. 7. The thermodynamic potential between the LOFF and 
normal state as a function of LOFF momentum. The thermo- 
dynamic potential is normalized to the magnitude of the ther- 
modynamic potential of the BCS state at 5 = Q. The momen- 
tum is normalized with respect to the chemical potential dif- 
ference. This curve corresponds to the value of the coupling 
l/kfO = -1.5 and 5^ = 0.72Zio- 



It is reassuring to note that the expression Eq.(57) for the 
thermodynamic potential is same as in Ref . [8] . This is also 
identical to the thermodynamic potential considered for 
quark matter in LOFF phase in Ref. [37] when degeneracy 
factors for the color and flavor are taken into account. 

Similar to Ref. [37], it can be shown that, for the case 
of a small gap as compared to the chemical potential mis- 
match parameter, i.e. z << 5, which will be the case near 
a second order phase transition to normal matter, the gap 
is given as [8,37] 



4(P? 



1)1-^ 



/ 3 = —{v - ^)'5/2 „ A [(j^ + 5)5/2 sf/A and, the momentum p = \q\/5^ is given as [8,37] 

15 4 15 V / 



(56) 



1 



-z\ 



Collecting all the terms, the thermodynamic potential Eq.(53jn the above, pc satisfies the equation 



IS given as 



27r2 2 



i-2 I « 

+ —V — 
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In- 
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+ {p-lfxl) 



(57) 



5a 



2^^ 
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Fig. 8. Super-fluid gap as a function of chemical potential 
difference. The coupling parameter here is taken as l/k/a = 
-1.5. 

for pc ~ 1.2. 

Here, Sc — 0.754zo is the maximum value of the chem- 
ical potential difference that can support the LOFF phase 
beyond which the system goes to the normal matter phase. 
For a general gap parameter which need not be small as 
compared to the chemical potential difference, one has to 
solve the gap equation Eq.(50) numerically. This is done 
for a given value of LOFF momentum p so that the ther- 
modynamic potential as given in Eq.(57) is calculated for 
a given value of the momentum. The minimization of the 
thermodynamic potential numerically over the magnitude 
of momentum gives the best value of the LOFF momen- 
tum. 

In the present calculations, however, we do not use the 
weak coupling gap equation, Eq.(50). Instead, for given 
values of coupling l/kfu, the chemical potential difference 
di, and LOFF momentum q, the coupled equations Eq.(23) 
and Eq.(26) are solved in a self consistent manner. Using 
the values of the gap and the average chemical potential 
so obtained we calculate the thermodynamic potential us- 
ing Eq.(26). For the numerical analysis, we first start with 
the weak coupling BCS regime Si, = Smax, the boundary 
that separates the LOFF phase and the normal phase. For 
Si, < Smax 1 we solve the gap Eq. (23) for a given value of 
|q|. Solution to the LOFF gap equations exist for a range 
of |q|. For each value of S we can determine the 'best 
|q|': the choice of |q| that has the lowest thermodynamic 



potential. In Fig. 7, the difference of thermodynamic po- 
tential between the LOFF phase and the normal matter 
is plotted as a function of the LOFF momentum for cou- 
pling l/kja — —1.5 and the chemical potential difference 
Si, = 0.72Z\o- We have normalized the thermodynamic 
potential with respect to the magnitude of the thermo- 
dynamic potential for the BCS phase. We note here that 
at this value of S, which is greater than the Clogston- 
Chandrasekhar limit, the normal matter has lower ther- 
modynamic potential than the homogeneous BCS phase. 
Thus the LOFF state is favored when the thermodynamic 
potential difference that is plotted in Fig. 7 is negative. For 
weak coupling this happens when Si, > 0.705Z\o, a value 
slightly lower than the Clogston-Chandrasekhar limit. In 
this region there are also other solutions of gap equation 
with homogeneous gap, but with higher value of the ther- 
modynamic potential. At each S < Smax we also compare 
the thermodynamic potential for the 'best |q|' and that 
of homogeneous BCS phase. We see that for S > Smin, 
LOFF state has lower thermodynamic potential as com- 
pared to the homogeneous BCS phase. At Smin there is 
a first order phase transition between the LOFF and the 
BCS phases. Thus in the 'LOFF window' <5max — ^mm, the 
LOFF state has lower energy as compared to the BCS or 
the normal phase. For weak coupling, our numerical values 
for the condensate, the momentum of the condensate and 
the LOFF window match with those obtained in Ref.[7, 
8]. 

The variation of the gap parameter as a function of 
difference in chemical potential is shown in Fig. 8 for cou- 
pling l/fc/a = — 1.5. Beyond Smax, the stress of difference 
in the chemical potentials leads to the vanishing of the 
condensate, as a second order phase transition. The varia- 
tion of the LOFF momentum with respect to the chemical 
potential difference is shown in Fig. 9. With y = |q|/fc/, 
(q, being the LOFF momentum), we have plotted here 
the ratio y/S as a function of the chemical potential dif- 
ference S which decreases monotonically from a value 1.29 
at S^/Ao = 0.704 to about 1.22 at S„/Ao = 0.754. These 
values are similar to those obtained in Ref.[8,37] in the 
weak coupling limit. 

Next we explore how the width of the LOFF win- 
dow depends upon the strength of the coupling. We see 
that from the weak coupling BCS side, as we increase the 
coupling the lower limit Smin of the LOFF window de- 
creases very slowly while the upper limit Smax increases. 
We might note here that the thermodynamic potential 
difference between the LOFF state and the normal mat- 
ter is very small. We had seen in the last subsection that 
the Clogston - Chandrasekhar limit remains almost con- 
stant in the BCS side of the Feshbach resonance (see e.g. 
Fig. 3). Smin which is very near to this value, thus is al- 
most constant on the BCS side of Feshbach resonance. 
With Smax increasing in this region, the LOFF window 
increases. However, the Clogston - Chandrasekhar limit 
increases with the coupling for positive scattering length 
and Smin of the LOFF window also increases resulting in 
reduction of the LOFF window in the positive coupling 
region. Thus the LOFF window turns out to be a non 
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Fig. 9. LOFF momentum in units of mismatch in chemical 
potential as a function of mismatch in chemical potential in 
units of gap at zero mismatch Aq. The coupling here is 1/fc/a = 
-1.5. 
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Fig. 10. The LOFF window in units of the gap at zero tem- 
perature as a function of coupling 1/fc/a 



monotonic function of the coupling and it becomes max- 
imum at 1/kfa ~ 0.04. Beyond this it decreases rapidly. 
For coupling greater than 1/kfa = 0.5; there is no longer 
any window of mismatch in chemical potential in which 
LOFF state occurs. However, these results can be modi- 
fied with a more involved ansatz of multiple plane waves 
as opposed to the single plane wave ansatz as considered 
in Eq.(3) or Eq.(7). Such multi plane wave ansatz was sug- 
gested to be favorable near the upper edge of the LOFF 
window [15]. In this context, it is worthwhile to mention 
also that the LOFF window can be considerably expanded 
if one considers optical lattices [35] rather than the pairing 
in free space as considered here. As discussed earlier, for 
couphngs 1/kfa > 1, the the interior gap state is preferred 
within a certain range of mismatch in chemical potential. 

4 Summary and discussions 

We have considered here a variational approach to dis- 
cuss the ground state structure for a system of fermionic 
atoms with mismatch in their fermi momenta. An explicit 
construct for the ground state is considered to describe 
two fermion condensate with nonzero momentum. The 
ansatz functions including the thermal distribution func- 
tions have been determined from the minimization of the 
thermodynamic potential. This is done for fixed chemical 
potentials. A similar analysis for fixed number densities 
leads to different results [33]. 

For comparison of thermodynamic potential, we have 
subtracted out the free energy of normal phase (with A ~ 
and q = 0) as in Eq.(lO). This difference is finite without 
introducing any arbitrary cut off in the momentum [14]. 
The four fermi coupling is also renormalized as in Eq.(19), 
by subtracting out the vacuum contribution and relating it 
to the s-wave scattering length as in Ref.s [10,24,31]. This 
makes all the quantities well defined and finite. Rewriting 
the gap equations and the number density equations for 
the symmetric case in terms of elliptic integrals as in Eqs. 
(29) and (30) is particularly useful regarding the numerical 
evaluation. 

We have not calculated here the Meissner masses [16, 
37] or the number susceptibihty [18] to discuss stability 
of different phases by ruling out regions of the parameter 
space for the gap function and the difference in chemical 
potentials. Instead, we have solved the gap equations and 
the number density equations self consistently and have 
compared the thermodynamic potentials. This apart, we 
have also extended the analysis to incorporate the possi- 
bility of condensate with finite momentum through the the 
ansatz for the ground state. In certain regions of chemical 
potential difference and couplings we have multiple solu- 
tions for the gap equation. In such cases we have chosen 
the solution which has the least value for the thermody- 
namic potential. 

In the weak coupling BCS limit we obtain the usual 
LOFF solution i.e. the condensate with finite momentum, 
in a small window of the chemical potential difference of 
about 0.05 times the gap at zero chemical potential dif- 
ference. This LOFF window increases as the coupling is 
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increased and becomes maximum at l/fc/a ^ 0.04 and 
vanishes at l/kfU = 0.5. Within the present approach we 
do not see any breached pairing solution being favored on 
the BCS side of the Fcshbach resonance. However, interior 
gap superfluid phase with a single fermi surface appears 
to be possible on the strong coupling BEC side of the 
Feshbach resonance with negative chemical potential sim- 
ilar to results in Ref.[18]. The transition between BCS to 
LOFF phase is a first order one with the order parameter 
varying discontinuously at Smin, while the LOFF phase to 
normal phase transition is second order at chemical po- 
tential difference S„iax- On the other hand, in the strong 
coupling BEC region, the phase transition from BCS to 
interior gap phase is a second order phase transition while 
the phase transition from interior gap phase to the nor- 
mal phase is a first order phase transition as the chemical 
potential difference^ is varied. 

These results are of course limited by the simplified 
ansatz as considered here giving a unified description for 
the uniform as well as spatially modulated superfluid. Sim- 
ilar ansatz with uniform condensates has been used earlier 
interpolating the BCS {g 0+, 1/a — oo) and BEC 
{g 00, 1/a 00 ) limits for the case of zero chem- 
ical potential difference [24,25]. The results as obtained 
here might nevertheless be regarded as a reference solu- 
tion with which other numerical or analytical solutions 
obtained from more involved ansatz for the ground state 
can be compared. 

Although the present numerical analysis has been done 
for zero temperature, the expressions to include tempera- 
ture effects have been derived here. Clearly, the effect of 
fluctuations involving the corrections from the collective 
modes will play an important role for high temperatures 
particularly for the strong coupling [24] . The effect of in- 
cluding different masses of the two species as well as using 
a realistic potential for the two atomic species will be inter- 
esting regarding the study of the phase structure. Further, 
the study of some of the transport properties in different 
phases including viscosity would be very interesting for 
cold atomic superfluids. Some of these calculations are in 
progress and will be reported elsewhere [38] . 
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